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Abstract 

Background: Water-soluble quinoprotein glucose dehydrogenase (PQQGDH-B) from 
Adnetobacter calcoacetiais has a great potential for application as a glucose sensor constituent 
Because this enzyme shows no activity in its monomeric form, correct quaternary structure is 
essential for the formation of active enzyme. We have previously reported on the increasing of the 
stability of PQQGDH-B by preventing the subunit dissociation. Previous studies were based on 
decreasing the entropy of quaternary structure dissociation but not on increasing the interaction 
between the two subunits. We therefore attempted to introduce a hydrophobic interaction in the 
dimeric internee to increase the stability of PQQGDH-B. 

Results: Amino acid residues Asn340 and Tyr4l8 face each other at the dimer interface of 
PQQGDH-B, however no interaction exists between their side chains. We simultaneously 
substituted Asn340 to Phe and Tyr4l8 to Phe or lie, to create the two mutants Asn340Phe/ 
Tyr4l8Phe and Asn340Phen'yr4l8lle. Furthermore, residues Leu280, Val282 and Val342 form a 
hydrophobic region that faces, on the other subunit, residues Thr4 1 6 and Thr4 1 7, again without 
any specific interaction. We simultaneously substituted Thr4 1 6 and Thr4 1 7 to Val, to create the 
mutant Thr4 1 6Val/Thr4 1 7Val. The temperatures resulting in lose of half of the initial activity of the 
constructed mutants were increased by 3-4''C higher over wild type. All mutants showed 2-fold 
higher thermal stability at 55°C than the wild-type enzyme, without decreasing their catalytic 
activities. From the 3D models of all the mutant enzymes, the predicted binding energies were 
found to be significantly greater that in the wild-type enzyme, consistent with the increases in 
thermal stabilities. 

Conclusions: We have achieved via site-directed mutagenesis the Improvement of the thermal 
stability of PQQGDH-B by Increasing the dimer Interface Interaction. Through rational design 
based on the quaternary structure of the enzyme, we selected residues located at the dimer 
Interface that do not contribute to the intersubunit interaction. By substituting these residues to 
hydrophobic ones, the thermal stability of PQQGDH-B was increased without decreasing its 
catalytic activity. 
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Background 

Water-soluble quinoprotein glucose dehydrogenase 
(PQQGDH-B) from Adnetobacter calcoaceticus has great 
potential for application as a constituent of an electron 
mediator-type glucose sensor. The conventionally utilized 
enzyme for glucose measurement, glucose oxidase 
(GOD), inherently utilizes oxygen as its electron acceptor 
during the oxidation of glucose. In contrast, PQQGDH-B 
is completely independent of oxygen, resulting in 
improved accuracy and rapidity during glucose measure- 
ment. However, because the substrate specificity and sta- 
bility of PQQGDH-B remain somewhat inferior to those 
of GOD, we have been engaged in the improvement of 
these enzymatic properties through protein engineering 
to further increase the application of PQQGDH-B in glu- 
cose monitoring systems [1-17]. 

The subunit structure of PQQGDH-B was determined to 
be homodimeric with no activity observed in its mono- 
meric form [17]. Correct quaternary structure is essential 
for the formation of active enzyme and dissociation of the 
dimer conformation triggers inactivation of the enzyme. 
We have previously reported on the increasing of the sta- 
bility of PQQGDH-B against the dissociation of quater- 
nary structure by chemical cross-linking [12], by 
constructing tethered enzyme [13], and by introducing 
Cys residues at the dimer interface to form a novel inter- 
subunit disulfide bond [14]. All of these attempts were 
based on decreasing the entropy by decreasing the possi- 
bility of dimer dissociation, but no attempts had been 
made to increase the interaction between the two 
subunits. 

In this paper, we report on our rational designing of 
hydrophobic interaction in the dimer interface to increase 
the stability of PQQGDH-B. We identified protein regions 
at the dimer interface where potential novel hydrophobic 
interactions could be introduced by amino acid substitu- 
tion, thereby increasing the stability. 

Results 

Modeling novel dimer hydrophobic core 

Among 19 amino acid residues located at the dimer inter- 
face, 8 residues were predicted not to be involved in the 
formation of hydrogen bonds, electrostatic interactions, 
or hydrophobic interactions (Fig. 1). We focused on the 
residues Asn340, Tyr418, Thr416, and Thr417, as they are 
not involved in the formation of the active site cavity and 
are therefore suitable candidates for amino acid substitu- 
tion. Although residues Asn340 and Tyr418 face each 
other on the surface of the dimer interface (Fig. 2), no 
interaction exists between their side chains. We therefore 
simultaneously substituted Asn340 to Phe and Tyr418 to 
Phe or He to create the mutants Asn340Phe/Tyr418Phe 
and Asn340Phe/Tyr418Ile, respectively. Furthermore, the 




Figure I 

The annino acid residues located at the dimer inter- 
face. The two subunits are represented in red and blue, 
respectively, using the RasMol molecular visualization soft- 
ware 26. The 19 amino acid residues at the interface are 
shown In space filling format, of which 8 residues (orange and 
light blue) are predicted not be involve in hydrogen bond for- 
mation, electrostatic interaction, or hydrophobic interaction 
at the interfiace. 



hydrophobic region composed of residues Leu280, 
Val282, and Val342 faces residues Thr416 and Thr417, 
again with no specific interaction. We therefore substi- 
tuted Thr4 1 6 to Val and Thr4 1 7 to Val to create the double 
mutant Thr416Val/Thr417Val. 

Characterization of mutant enzymes 

The activity and stability of each of the three constructed 
double mutant enzymes were then analyzed. All three 
mutant enzymes showed slightly higher thermal stability 
than wild-type PQQGDH-B upon incubation for 10 min 
at various temperatures (Fig. 3). The temperatures result- 
ing in lose of half of the initial activity were shifted by 
approximately 3°C higher in the mutants compared to 
the wild type (Table 1). As the time course of thermal inac- 
tivation at SS^C follows first-order kinetics (Fig. 4), half- 
lives were calculated using logarithmic regression of resid- 
ual activity. The wild-type enzyme inactivates at 55 °C 
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Table I : Kinetic parameters and thermal stability of PQQGDH-Bs. 





/max(U/mg) 


fCm(mM) 




t,a(min)'' 


Wild type 


3030 


20 


53.9 ± 0.9 


9.5 ± 1.4 


Asn340PheATyr4l8Phe 


3100 


20 


57.7 ± 0.4 


14.9 ± I.I 


Asn340Phe/r/r4l8lle 


2500 


20 


57.5 ± 0.3 


15.5 ± 1.5 


Thr4l6Vain-hr4l7Val 


2800 


16 


56.5 ± 0.9 


14.8 ± 1.4 



e at which half the initial activity is lost in 1 0 min. 
12 represents the half-life at 55°C. 




Figure 2 

Hydrophobicity of PQQGDH-B din 

Hydrophobic regions are shown in green and hydrophilic 
regions are shown in blue. Residues that have been substi- 
tuted in this study are indicated. The structural images were 
generated using Molecular Operating Environment 




20 25 30 35 40 45 50 55 60 65 70 75 80 8 
10-minute treatment temperature (°C) 

Figure 3 

Thermal stability of wild-type and mutant 
PQQGDH-Bs. Residual activity was measured at 25°C after 
lO-min incubations at different temperatures of the following 
protein samples (0.075 \i^mL): Wild-type ♦, Asn340Phe/ 
Tyr4l8Phe □, Asn340Phen'yr4 1 8lle A, and Thr4l6Val/ 
Thr4l7ValO. 



with a half-life of 9.5 + 1.4 min, while all three double 
mutants showed greater thermal stability, with half-lives 
of 5-6 °C greater (Table 1). 

Investigation of the kinetic parameters of the wild-type 
and mutant enzymes shows that the mutations did not 
significantly affect the overall kinetic properties of 
PQQGDH-B (Table 1). The specific activities of 
Asn340Phe/Tyr418Phe (3100 U/mg), Asn340Phe/ 
Tyr418Ile (2500 U/mg), andThr416Val/Thr417Val (2800 



U/mg) were close to that of the wild-type enzyme (3030 
U/mg). Except forThr416Val/Thr417Val, which had a Km 
value of 16 mM, the mutants had Km values identical to 
20 mM Km value of the wild-type enzyme. 

Discussion 

PQQGDH-B is a 6-blade P-propeller protein, with each 
blade consisting of a 4-stranded anti-parallel P-sheet (W- 
motif) [18]. The strands in each W-motif are labeled A-D 
from the inside to the outside of the molecule [18,19]. All 
strands are connected by loops, which are named accord- 
ing to the strands they connect. Based on PQQGDH-B 
structural information, these loop regions have been 
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Figure 4 

Time course of thermal inactlvation of wild-type and 
mutant PQQGDH-Bs. Protein samples (0.075 p.g/mL) 
were incubated at SS'C and aliquots were taken at different 
times to measure residual activit/. The analyzed samples con- 
tained the following PQQGDH-Bs: Wild-type ♦. Asn340Phe/ 
Tyr4 1 8Phe □, Asn340Phe/Tyr4 1 8lle A, Thr4 1 6Val/ 
Thr4l7ValO. 



Table 2: Predicted binding en 


ergy of each mutan 


t subunit 


calculated using AMBER89 a 


d CHARMM22 fbro 


e fields. 




Binding energy (l(cal/mol) 




AMBER89 


CHARMM22 


Wild type 


-249 


-118 


Asn340PheATyr4l8Phe 


-291 


-134 


Asn340PheAryr4l8lle 


-273 


-126 


Thr4l6Val/Thr4l7Val 


-278 


-131 



associated with a number of important functions, such as 
substrate binding, co-factor binding, and formation of the 
enzyme active site. As with other P-propeller proteins, the 
catalytic site and substrate-binding pocket of PQQGDH-B 
is made up of the cleft formed by loops DA and BC [20]. 
The enzyme surface composed of loops AB and CD is 
therefore located opposite the functional region. In the 
present study, we have introduced mutations at Asn340, 
Tyr418, Thr416, and Thr417, which are all located in the 
loop 5CD region. Considering that these residues are 
located far from the functional region and do not contrib- 



ute in the structure of functional region, it is not surpris- 
ing that their substitutions did not significantly alter the 
enzyme's kinetic parameters, particularly its catalytic 
activity. 

Based on the predicted structure of the mutant enzymes 
shown in Figure 5, the binding energy of each subunit was 
calculated, using both the AMBER89 and CHARMM22 
force fields, and compared with those of the wild-type 
enzyme (Table 2). As expected, the increases in hydropho- 
bicity at the subunit interface resulted in increases in their 
binding energies calculated by both methods. Estimation 
of the number of hydrophobic interactions based on these 
same predicted models revealed 4 to 5 novel hydrophobic 
interactions at the interface of all the mutants, while none 
were found in the wild-type one. These results based on 
structural predictions are consistent with the observed 
improvements in thermal stability. 

The addition of hydrophobic interactions in other 
enzymes has been reported to result in thermostability 
increases of 2 to 10°C [21-24], comparable to the results 
of the current study. Recent observations of enzymes from 
thermophilic organisms indicate that their extraordinary 
thermal stability is due to hydrophobic interactions. Oli- 
gomeric enzyme stability was also reported to be 
improved through engineering to increase the hydropho- 
bic interaction at the oligomer interface [25]. However, 
these studies were based on homology analyses between 
mesophilic and thermophilic bacteria together with ran- 
dom mutagenesis library analyses [21-24]. Although 
some sequences have been found to be homologous to 
PQQGDH-B, these are all putative ORFs with no func- 
tional information reported. Therefore, there exists no 
reliable template to help improve the stability of this 
enzyme by increasing the dimer interface interaction. 

Conclusions 

We improved PQQGDH-B's thermal stability by increas- 
ing the dimer interface interaction through rational design 
based on its quaternary structure. We demonstrated that 
this can be achieved by selecting residues located at the 
dimer interface that do not contribute to the intersubunit 
interaction and substituting them to hydrophobic ones 
via site-directed mutagenesis. In each case tested, the 
enzyme's thermal stability was increased without decreas- 
ing its catalytic activity. This rational design approach will 
provide relevant information for future designs by com- 
bining with other mutant PQQGDH-Bs with narrowed 
substrate specificity and improved catalytic efficiency. 

Methods 

Chemicals 

Glucose, phenazine methosulfate (PMS), 2,6-dichloroph- 
enolindophenol (DCIP), and magnesium chloride were 
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Figure 5 

Hydrophobicity of mutant PQQGDH-B dimer interface. Hydrophobic regions are shown in green whiie hydrophiiic 
regions are shown In blue. The interfaces shown are those of the Asn340Phe/Tyr4 i 8Phe (ieft), Asn340Phe/Tyr4i8iie (middie), 
and Thr4l6Val/Thr4l7Val (right) mutants of PQQGDH-B, with their respective mutation sites circied. The structural images 
were generated using Moiecuiar Operating Environment 



obtained from Kanto Kagaku (Tokyo, Japan), 3-(N-mor- 
pholino) propane sulfonate (MOPS) from Dojin 
(Kumamoto, Japan), and pyrroloquinoline quinone from 
Mitsubishi Gas Chemical Company (Tokyo, Japan). All 
other regents were of analytical grade. Kpnl was obtained 
from TOYOBO (Osaka, Japan) and Hindlll from New 
England BioUbs (Beverly, USA). 

Site-directed mutagenesis 

The structural gene for wild-tj^e PQQGDH-B was previ- 
ously amplified by polymerase chain reaction (PGR) and 
inserted into the expression vector pTrc99A (Pharmacia) 
to create pGB [15]. A 1.2-kbp KpnI-Hindlll fragment con- 
taining the PQQGDH-B gene was transferred firom pGB to 
pKFlSk and mutagenesis was carried out with the Mutan- 
Express Km kit (Takara) according to the manufacturer's 
instructions with the oligonucleotides Asn340Phe (5- 
GGTGGGACAAAGMTrrACCAGTCC-3'), Tyr418Phe (5'- 
CGGTACAGCGTCATCAAAAGTAGTGC-3'), Tyr418Ile 
(5 '-GGGTACAGCGTCATCMrAGTAGTGC-3 '), and 
Thr416Val/Thr417Val (5'-CACCGTCATCAT- 
A AACAAC GCTATAAGTrGGATC-3'). The mutations 
(underlined) were confirmed by automated DNA 
sequencing (ABI PRISM Genetic analyzer 310, Applied 
Biosystems). The mutated genes were digested with Kpnl 
and Hindlll and were replaced into pGB to construct 
expression vectors containing mutated PQQGDH-B. 
Numbering of the amino acid positions starts from the 
first residue of the signal peptide (24 residues). 



Enzyme preparation and assay 

The PQQGDH-B genes were expressed in Escherichia coli 
and the enzymes purified as previously reported [15,16]. 
Following a 30-min preincubation in 10 mM MOPS- 
NaOH (pH 7.0) containing 1 ^M PQQ and 1 mM CaCl^ 
at room temperature (25 °C) to produce the holoenzyme, 
GDH activity was measured by using 0.6 mM PMS and 
0.06 mM DCIP. The enzyme activity was determined by 
measuring the decrease in absorbance of DCIP at 600 nm. 

Analysis of PQQGDH-B stability 

The thermal stability of wild-type and mutant PQQGDH- 
B was determined with 0.075 (ig/mL protein, as previ- 
ously reported [15]. Thermal inactivation experiments 
were carried out by incubating each holoenzyme in 200 
jiL of 10 mM MOPS-NaOH, pH 7.0, at 55 "C. Aliquots 
were taken every 5 min and kept at 4°C for 2 min, 
followed by incubation at room temperature for 30 min. 
The residual enz3mie activity was determined as described 
above. Since the initial time course for thermal inactiva- 
tion at 55''C followed first-order kinetics, the thermal sta- 
bility of each mutant enzyme was expressed as a half-life. 
The thermal stability of Asn340Phe/Tyr418Phe, 
Asn340Phe/Tyr418Ile and Thr416Val/Thr417Val were 
also determined by incubating purified enzyme at various 
temperatures for 10 min. Hie residual activities were 
determined as described above, and were compared with 
the initial activities. 
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Prediction of three-dimensional structure and quaternary- 
dimensional structures 

Three-dimensional and quaternary structures were pre- 
dicted using Molecular Operating Environment (MOE) 
(Chemical Computing Group Inc., Quebec, Canada). By 
using the available PDB data of the wild-type PQQGDH- 
B, IQBI [18], we made the appropriate substitutions with 
all possible side-chain orientations to predict the struc- 
tures of the Asn340Phe/Tyr4 1 8Phe, Asn340Phe/Tyr4 1 Bile 
and Thr416Val/Thr417Val mutants. After addition of 
hydrogen atoms to the PQQGDH-B structure and optimi- 
zation of orientation of some hydrogen atoms by MOE, 
the structures were subjected to energy minimization 
using the AMBER89 or CHARMM22 force field within the 
MOE program until the final energy gradient was < 0.01 
kcal/mol- A. 
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